Heavy-quark transport coefficients in a hot viscous quark-gluon plasma medium 
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The heavy-quark (HQ) transport coefficients have been estimated for a viscous quark-gluon plasma 
medium, utilizing a recently proposed quasi-particle description based on realistic QGP equation of 
state (EoS). Interactions entering through the equation of state significantly suppress the tempera- 
ture dependence of the drag coefficient of QGP as compared to that of an ideal system of quarks and 
gluons. Inclusion of shear and bulk viscosities through the corrections to the thermal phase space 
factors of the bath particles alters the magnitude of the drag coefficient and the enhancement is sig- 
nificant at lower temperatures. The competition between the effects of realistic EoS and dissipative 
corrections through phase space factor the former eventually dictate how the drag coefficient would 
behave as a function of temperature, and how much quantitatively digress from the ideal case. The 
observations suggest significant impact of both the realistic equation of state, and the viscosities, 
on the HQs transport at RHIC and the LHC collision energies. 
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I. INTRODUCTION 

Experimental heavy-ion collision (HIC) program at 
Relativistic Heavy Ion Collider (RHIC) indicates the pro- 
duction of a liquid like state of the matter where the 
properties of the system is governed by quarks and glu- 
ons. Such a state of matter is called as quark-gluon 
plasma (QGP). Two of the most striking finding of the 
RHIC which lead to the imposition of strongly interact- 
ing liquid-like behaviour are viz., the large elliptic flow 
shown by QGP, and the jet quenching phenomenon [l[. 
Preliminary results from LHC HHII confirm the above 
mentioned observations regarding the QGP. The strongly 
coupled picture of the QGP is found to be consistent with 
the lattice simulations of the hot QCD equation of state 
(EoS) which predict a strongly interacting behavior 

even at temperatures which are of the order of a ~ 2T C , 
where T c is the temperature for the quark-hadron tran- 
sition. Therefore, it calls for the adaption of a realistic 
equation of states for the QGP (such as obtained form 
lattice QCD) in the theoretical and phenomcnological in- 
vestigations regarding its bulk and transport properties. 
We shall closely follow this point of view in the present 
analysis in the later sections. 

On the other hand, hadrons containing HQs (c, c and 
b, b) are of great interest in investigating the properties of 
the QGP, since their physical properties get significantly 
modified while traveling through the QGP medium. This 
fact has been reflected in the particle spectra at RHIC, 
and the LHC energies. Further, HQ thermalization time 
is larger than gluons and light quarks, and they do not 
constitute the bulk of the QGP. Since their formation 
occurs in the early stages of the collisions, they can travel 
through the thermalized QGP medium, and can retain 
the information about the interaction with them very 



effectively. For instance, it is pertinent to ask whether a 
single cc can stay together long enough to form a bound 
state (say J/tp) at the hadronization state. To address 
this, one requires to describe the dynamics of the HQs 
propagating through the QGP. 

Therefore, one can envisage the HQ transport in the 
QGP medium as follows. The non-equilibrated HQs can 
travel in the equilibrated QGP medium, and one has to 
deal theproblem within the framework of Langevin dy- 
namics [7[ • This is to say that the HQs undertake random 
motion in the heat bath of the equilibrated QGP. Recall, 
that the QGP goes through a hydrodynamic evolution 
process before it reaches the hadronization and subse- 
quently the hadrons freeze-out. The pertinent question 
to ask is, whether a HQs maintain equilibrium during 
this entire process of evolution or not. It has been ob- 
served @] within the ambit of Langevin dynamics and 
pQCD that they do not maintain the equilibrium in the 
RHIC and the LHC energies. 

The most important observable which encode the 
medium effects carried with them by the HQs while trav- 
eling in the QGP, is the nuclear modification factor, 
Raa- It has been observed that their energy loss in the 
QGP due to gluon radiation is insufficient to describe the 
medium modification of the spectrum ■ Therefore, 

one has to look at the collisions since they have differ- 
ent fluctuation spectrum than radiation, and might con- 
tribute significantly as one thought off initially 12[ . The 
collisional effects can be captured well in the HQ drag and 
diffusion coefficients which have been calculated within 
weak coupling QCD by several authors. The formalism, 
and details are offered in Sec. IIA. 

The temperature, T and chemical potential, (j,b depen- 
dence of the drag and diffusion coefficient enter through 
the thermal distributions of gluons and light quarks. In 
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the present case, we ignore the [1b dependence in view 
of the fact that the QGP produced at RHIC and LHC 
energies at the mid-rapidity region has negligibly small 
net baryon density. Therefore, one has to implement the 
realistic QGP EoS in terms of appropriate form of the 
thermal distribution functions. Another important as- 
pect is the viscosities (both shear and bulk) of the QGP. 
One also has to include the viscous modifications in the 
distribution functions by maintaining the near equilib- 
rium picture so that the mathematical formalism of HQ 
transport remain intact. This is possible since the QGP 
posses small shear and bulk viscosities for temperature, 
T > T c (T c is the quark-hadron transition). Both these 
aspect have been included for the first time while study- 
ing the temperature dependence of HQ drag and diffusion 
coefficients in the QGP medium. Discussion regarding 
the former is offered in Sec. IIB, and the latter in IIC. 
Our results suggest that both of these modifications have 
significant impact on the temperature dependence of HQ 
transport coefficients. Interestingly, the inclusion of in- 
teraction through the EoS modulate the transport coef- 
ficients quite significantly as compared to those obtained 
by employing an ideal EoS. 

The paper is organized as follows. In Sec. IIA, formal- 
ism of HQ drag and diffusion has been presented. Sec. 
IIB deals with the quasi-particlc description of hot QCD 
and modeling of thermal distribution functions of gluons, 
and quarks which constitute hot QGP medium. Sec. IIC 
deals with the viscous modifications to the thermal distri- 
bution functions. In Sec III, results and discussions have 
been presented for the drag and diffusion coefficients with 
the realistic EoS. In Sec. IV, this has been investigated 
including the viscous modifications. Sec. V is devoted to 
summary and conclusions. 



II. HEAVY-QUARK DRAG AND DIFFUSION 
IN THE VISCOUS QGP 



A. Formalism of HQ drag and diffusion 

The special role of HQs to characterized QGP stem 
from the fact that they are produced very early in the 
collisions and remain extant throughout the evolution 
and hence can witness the entire evolution of the sys- 
tem. Moreover, HQs do not constitute the bulk part of 
the system because their mass is considerably larger than 
the temperature attainable in HIC at RHIC and LHC. 

In the present work, we consider the elastic interac- 
tion experienced by HQs when propagating in the QGP 
medium. For the process, c(p) + l(q) — > c(p') + l(q') (I 
stands for gluon and light quarks and anti quarks), the 
drag 7 can be calculated by using the following expres- 
sion (SHI: 



where A4 is given by 
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gq being the statistical degeneracy of the HQ propagat- 
ing in the medium. The above expression indicates that 
the drag coefficient is the measure of the thermal aver- 
age of the momentum transfer, p — p' due to interaction 
weighted by the square of the invariant amplitude, | M \ 2 . 
The factor f(q) denotes the thermal phase space for the 
particle in the medium. We shall see in the subsequent 
subsection that /(g) will involve three types of thermal 
phase space distribution corresponding to the gluons, 
light-quarks (up and down) and the strange quarks and 
their anti-quarks. Therefore, f(q) jointly denote these 
three sectors as, 



/(<?) = {fgJgJs} 



(3) 



in the absence of the dissipative effects coming through 
the shear and the bulk viscosity. It may be recalled here 
that for zero baryonic chemical potential the thermal 
phase space for quarks and anti-quarks are same. In the 
presence of dissipation, we denote it as, 



fin) — {fgni fqqi fss} 



(4) 



We shall discuss these distributions in the subsequent 
subsections. 

Similar to drag the diffusion coefficient, Bq can be de- 
fined as: 
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With an appropriate choice of Tip') both the drag and 
diffusion co-efficients can be evaluated from the following 
expression: 
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7 = PiAi/p 2 



(1) 



wheres is the Mandelstam variable, Eq = ^Jp 2 + jtiq is 

the energy of the HQ, E q is the energy of the bath particle 
w(a,b,c) = a 2 + b 2 + c 2 — 2ab — 2bc — 2ac, is the triangular 
function. 

To distinguish among various cases, following nota- 
tions will be utilized in the subsequent sections. We shall 
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use the notations, 7^ and 7^ for the the drag coefficient 
of the HQs in QGP background with ideal EoS in the 
presence of shear and the bulk viscosities respectively. 
On the other hand in the realistic QGP EoS, the same 
quantities will be denoted as 7 eos and r y eos respectively. 
The ideal counter part (no viscous effects and ideal EoS 
for QGP) will be denoted as jid- Whenever both 77 and 
C are taken in to account, the drag coefficients will be 
denoted as 7^+^ in the case of ideal QGP background, 
and as 7^+ S £ in the case of realistic QGP EoS. 



v q = 2 x 2 x N c x 2 for light quarks, v a = 2 x 2 x N c x 1 for 
the strange quark for SU(N C ). Here, we are dealing with 
SU(3), so N c — 3. Since the model is valid in the de- 
confined phase of QCD (beyond T c ), therefore, the mass 
of the light quarks can be neglected as compared to the 
temperature. 

Clearly, the determination of f g , f q , and f s require 
the temperature dependence of the free parameters z g 
and z q respectively. To do this, we rewrite lattice QCD 
pressure symbolically as, 



B. The quasi-particle description of hot QCD 



P = P a +Pc 



(8) 



Now we discuss the quasi-particle nature of the hot 
QCD medium which is employed recently in [iTi . This 
description has been developed in the context of the 
recent (2+l)-lattice QCD EoS [l5| at physical quark 
masses. There are more recent lattice results with the 
improved actions and refined lattices 0, for which we 
need to re-look the model with specific set of lattice data 
specially to define the effective gluonic degrees of free- 
dom. This is beyond the scope of the present analysis. 
Henceforth, we will stick with the one set of lattice data 
which is utilized in the model [l6[. 

The model initiates with an ansatz that the Lattice 
QCD EoS can be interpreted in terms of non-interacting 
quasi-partons having gluon and quark effective fugaci- 
ties z g and z q respectively which encode all the interac- 
tion effects. In this approach, the hot QCD medium is 
divided into two sectors, viz., the effective gluonic sec- 
tor, and the matter sectors (light and strange quark). 
The former refers to the contribution of gluonic action 
to the pressure which also involves contributions from 
the internal fcrmion lines. On the other hand, latter in- 
volve interactions among quark, anti-quarks, as well as 
their interactions with gluons. The ansatz can be trans- 
lated to the form of the equilibrium distribution func- 
tions, f eq = {fg,fq,f s } (this notation will be useful later 
while writing the transport equation in both the sector 
in compact notations) as follows, 



h = 



Is = 



Zg exp(-/3p) 
1 - z g exp(-/3p) 

Zgexp(-/3p) 
1 + z q exp(-/3p) 



z q exp(-f3y/p 2 + m 2 ) 



1 + z q exp(— (3y/p 2 + m 2 ) 



where p = \p\, m s denotes the mass of the strange 
quark(which we choose to be O.lGeV), and (3 = T^ 1 
denotes inverse of the temperature. We use the nota- 
- 1) for gluonic degrees of freedom, 



where P is full (2+l)-flavor lattice QCD pressure, P g is 
the contribution coming from the gluonic action alone 
and Pqg is the remaining part after subtracting P g from 
P. The fraction P g , is utilized to define effective gluonic 
(quasi-gluons) degrees of freedom, and P qg to effective 
quark-antiquark degrees of freedom through an effective 
Grand canonical partition function Z 
relation P/3V = ln(Z) as follows, 



Zg x Zqg via the 



^ J d 3 pln(l-z s cxp(~/3p)) 
= (P - P g )fiV = £| / d 3 p\n(l + ZqcM-Pp)) 



p 9 py 

+ S J d 3 pln(l + Zqexp(-(3^/p 2 + m 2 )). (9) 

One can read-off Z g and Z qg which arc their standard 
definitions for bosons and fermions. We know the left 
hand side of these equation from lattice QCD, and by 
calculating roots of these equations at each temperature 
in lattice QCD, we can determine z g and z q as a function 
of temperature. This has been done exactly numerically. 
Note that we are working at zero baryon chemical poten- 
tial, therefore quark and antiquark distribution function 
are the same. This fact has been taken care of by count- 
ing the anti-quark degrees of freedom in v q and v s . 

It is worth emphasizing that the effective fugacity is 
not merely a temperature dependent parameter which 
encodes the hot QCD medium effects. Its physical signifi- 
cance is reflected through the modified dispersion relation 
both in the gluonic and matter sector obtained from the 
thermodynamic relation of energy density e = — dp ln(Z). 
One thus find that the effective fugacities modify the sin- 
gle quasi-parton energy as follows, 



uj g = p + T 2 dTln(z g ) 
Uq — p + T 2 dTln(z q ) 
uj s = \/p 2 + m 2 + T 2 dTln(z q ), 



(10) 



tion Vg = 2 (AT 



and this lead to the new energy dispersions for gluons 
(u! g ), light-quark antiquarks (u) q ), and strange quark- 
antiquarks, (u) s ). These dispersion relations can be expli- 
cated as follows. The single quasi-parton energy not only 
depends upon its momentum but also gets contribution 
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from the collective excitations of the quasi-partons. The 
second term is like the gap in the energy-spectrum due 
to the presence of quasi-particlc excitations. This makes 
the model more in the spirit of the Landau's theory of 
Fermi -liquids. A detailed discussions of the interpreta- 
tion and physical significance of z g , and z q , we refer the 
reader to UM . 



C. Viscous QGP: modification to thermal 
distributions 



equally significant as the shear viscosity for the temper- 
atures which are close to T c [22l - [25| . This is attributed 
to the large value of the interaction measure, and its be- 
havior with the temperature near T c 122)' . 

In the present analysis the mass of the strange quarks 
does not play any prominent role at high temperature 
(T > T c ) because m/T < 1. In fact, the mass effects are 
negligibly small (of the order of (m/T) 2 . Therefore, for 
the ease in the numerical evaluation the strange quark 
mass can be ignored. 



Viscosities usually modifies the equilibrium thermal 
distributions as follows, 



f = fe q +5f n +5f Q . 



(11) 



Here, / jointly denote the viscous modified thermal dis- 
tribution for gluons and quark-antiquarks (in the present 
case gluons, light quarks and strange quark). Therefore, 
/ = {fgn, fqm fqs}- For the viscous corrections 6f v , Sf c 
we utilize the results of [171 ] by using only first order terms 
in the expansion of shear and bulk part of the stress ten- 
sor, and choosing the local rest frame of the fluid. Pre- 
cisely, we have minimally extended the expressions in the 
present case by substituting the thermal distributions of 
quarks and gluons as the form obtained in the effective fu- 
gacity quasi-particle model. This is solely based on that 
the quasi-partons in the present case are non-interacting 
up to their respective effective fugacities that encode all 
the interactions. This leads to the following expressions, 



fgn — fg + 



fg(l + fg) \, P 2 



9 ' 2T 3 T 



2\ r l 
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fq^-fq) \,V 2 



i ' 2T 3 t 



5 s 

fs + 



2T 3 1 



■P,)- 



P Z )- 



5 s 



(12) 



(13) 



(14) 
(15) 



where r\ and C are the shear and bulk viscosities of the 
QGP and s is the entropy density. We will employ the 
same values of rj/s and £/s while comparing the effects of 
ideal and non-ideal EoS with viscosities. For simplicity, 
we assume same values for rj/s and (/s for gluons, and 
quarks. Here, r is taken to be the thermalization of the 
QGP. 

There are several interesting reports on the estima- 
tion of a small value of rj/s for QGP, and the emergence 



of near perfect fluid picture [18H2l|. It has been real 



III. RESULTS AND DISCUSSIONS 

Here, we show the temperature dependence of the HQ 
drag coefficient employing the realistic EOS for the QGP 
and the dissipative environment mainly induced by the 
shear and bulk viscosities. Our focus is mainly on the the 
impact of the the viscous effects, and the realistic QGP 
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ized in the recent past that bulk viscosity of the QGP is 



FIG. 1: 7/d stands for the drag coefficient for ideal case and 7^ 
stands for drag coefficients when the effect of shear viscosity 
is included through the phase space of the bath particles. The 
variation of the ratio of drag coefficients as a function of tem- 
perature with different values of rj/s at a given thermalization 
time, r are shown. 



Our aim is served by the ratio in the presence of vis- 
cosities and the EOS with respect to that is obtained by 
assuming the ideal EOS for the QGP medium in which 
the heavy quarks propagate. These effects arc shown in 
Figs. EE 

In Fig [T] the effect of shear viscosity on drag coefficient 
has been depicted as a function of temperature through 
the ratio between the drag coefficients with shear viscos- 
ity to drag coefficient with ideal case. The drag coefficient 
for zero shear viscosity (ideal case) is higher compared to 
the case of nonzero shear viscosity in the entire temper- 
ature range considered here. Drag coefficient decrease 
with increase of rj/s. The same quantity is displayed in 
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FIG. 2: Same as Fig. [JJfor i)/s — l/(47r) at two values of r. FIG. 4: Same as Fig. 0for £/s = 0.01 at two different values 
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FIG. 3: 7^ stands for drag coefficients when the effect of bulk 
viscosity is included through the phase space of the bath par- 
ticles. The variation of the ratio of drag coefficients as a 
function of temperature with different values of (/s at a given 
value of thermalization time r are depicted. 



2 
1.8 
1.6 
>- 1.4 

+ 



1.2 - 
1 - 
0.8 - 



r|/s=0.08, £/s=0.02, x=0.3 

-- r|/s=0.08, £/s=0.01,x=0.3 
n/s=0.16, C/s=0.01,x=0.3 




0.2 



0.4 

T (GeV) 



0.6 



0.8 



FIG. 5: fri+c stands for drag coefficients when the effects 
of both the shear and bulk viscosities are included through 
the thermal distribution of the bath particles in evaluating 
the drag coefficient. The variation of the ratio of the drag 
coefficients as a function of temperature with different values 
of C,/s,rj/s and r are shown. 



Fig [5] as a function of temperature for different values of 
thermalization time at a given value of rj/s. It is observed 
that effect of shear viscosity will decrease if we increase 
the thermalization time. In Fig [3] the effect of bulk vis- 
cosity on drag coefficient has been studied as a function 
of temperature taking a ratio between the drag coeffi- 
cients with bulk viscosity to drag coefficient with ideal 
case. It is found that bulk viscosity increases the drag 
coefficient. If we increase the (/s, drag coefficient is also 
enhanced as shown in Fig [3] In Fig [4] the effect of ther- 
malization time is displayed. It is found that the effect of 
bulk viscosity will decrease if we increase the thermaliza- 
tion time. Taking the effect of both the shear and bulk 
viscosity into consideration the drag coefficient has been 
calculated (7 r) +c) and its variation with the temperature 



is displayed in Fig [5] As drag coefficients decrease with 
i)i ' s and increase with £/s, it is a interplay between them 
when both these viscosity effects are included simulta- 
neously. It is found that the low temperature domain 
is dominated by bulk viscosity where as high temper- 
ature part is dominated by shear viscosity and that is 
why at low temperature drag coefficient is increased and 
high temperature drag coefficient is reduced as shown in 
Fig [5] Effect of the EoS on drag coefficient is plotted in 
Fig [5] by taking a ratio with the ideal case with respect 
to temperature. It is observed that effect of the EoS sup- 
press the drag coefficient at low temperature. In Fig [7] 
effects of both shear and bulk viscosities along with the 
EoS have been presented as a function of temperature at 
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FIG. 6: 7 eos stands for drag coefficients when effect of EoS is 
considered in evaluating the drag coefficient. The variation of 
the ratio of the drag coefficients as a function of temperature 
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FIG. 7: stands for drag coefficients when the effects of 

both the shear and bulk viscosity as well as the effect of EoS 
are included in evaluating the drag coefficient. The variation 
of the ratio of the drag coefficients as a function of tempera- 
ture is depicted here. 



given values of r, rj/s, £/s. The reduction in the drag is 
significant at the temperature domain relevant of nuclear 
collisions at RHIC and LHC energies. Therefore, this re- 
sults will have crucial effects on the nuclear suppression 
of flavors measured at these energies. 

The effects of phase space corrections due to non-zero 
shear and bulk viscosities and realistic EoS on diffusion 
coefficient (Bo) has not been shown explicitly here. How- 
ever, the effects can be estimated from the drag coeffi- 
cient (7) by using the Einstein relation: Bq = jmT. 
IV. SUMMARY AND CONCLUSIONS 



The drag and diffusion coefficients of the HQs propa- 
gating through a thcrmalizcd system of quarks and glu- 
ons have been evaluated by incorporating the corrections 
to the phase space factors arising from the non-zero vis- 
cosities of the system which restrain the system slightly 
away from the equilibrium. This corrections seem to be 
significant for the temperature range expected to be at- 
tained in the system produced at HIC at RHIC and LHC 
collision energies. The effect of realistic EoS on the drag 
has also been ascertained. The combined effects of vis- 
cosities and the EoS reduces the drag significantly. These 
results will have crucial consequences on the heavy ion 
phenomenology at RHIC and LHC energies. 
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